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Abstract 


Hydrological dimensioning or revision of the hydraulic works and the 
elaboration of flood risk maps is based on the so-called Design Floods, 
which are maximum flows of the river associated with low probabilities of 
exceedance. The most reliable way of estimating such predictions is 


through Flood Frequency Analysis (FFA). Its fundamental assumption is 


Tecnología y ciencias del agua, ISSN 2007-2422, 12(3), 164-203. DOI: 10.24850/j-tyca-2021-03-05 


2021, Instituto Mexicano de Tecnología del Agua 


Tecnología y Open Access bajo la licencia CC BY-NC-SA 4.0 
Ciencias¿Agua (https://creativecommons.org/licenses/by-nc-sa/4.0/) 


that the stochastic process under study is stationary, ¡.e., it does not 
change withtime. Construction of reservoirs, urbanization and changes in 
land use in the basin, as well as global or regional climate change, alter 
the hydrological processes and generate records of annual flows that are 
non-stationary, showing trends and changes in variability. For the FFA of 
such registries, the extreme value theory has been extended to apply ¡ts 
classical distribution, the Generalized Extreme Values (GEV), with 
parameters of location (u) and scale (a) varying with time (t), which is 
introduced as a covariate. In this work the L-moment method is 
presented for the fit of the probabilistic model GVE11 whose parameters u 
and a vary linearly over time. Three numerical applications are described. 
Conclusions highlight the simplicity of the exposed method and its 
importance in the estimation ofthe expected predictionsin non-stationary 


annual maximum data senies. 


Keywords: Covariate, L-moments, GEV distribution, standard error of 


fit, trend, linear regression, residuals. 


Resumen 


El dimensionamiento o la revisión hidrológica de las obras hidráulicas y la 
elaboración de mapas de riesgo por inundación se realizan con base en 
las llamadas crecientes de diseño, que son gastos máximos del río 
asociados con bajas probabilidades de excedencia. La manera más 
confiable de estimar tales predicciones es a través del Análisis de 
Frecuencias de Crecientes (AFC), cuya suposición fundamental es que el 


proceso estocástico bajo estudio es estacionario, es decir, que no cambia 
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con el tiempo. La construcción de embalses pequeños, la urbanización y 
los cambios de uso del suelo en la cuenca, así como el cambio climático 
global o regional alteran los procesos hidrológicos y generan registros de 
crecientes anuales que son no estacionarios, al mostrar tendencias y 
cambio en su variabilidad. Para el AFC de tales registros se ha extendido 
la teoría de valores extremos, a fin de aplicar su distribución clásica, la 
General de Valores Extremos (GVE), con parámetros de ubicación (u) y 
escala (a) variables con el tiempo (t), que se introduce como covariable. 
En este trabajo se expone el método de momentos L para el ajuste del 
modelo probabilístico GVE11, cuyos parámetros u y a varían linealmente 
con el tiempo. Se describen tres aplicaciones numéricas. Las conclusiones 
destacan la sencillez del método expuesto y su importancia en la 
estimación de las predicciones buscadas en series de datos máximos 


anuales no estacionarios. 


Palabras clave: covariable, momentos L, distribución GVE, error 
estándar de ajuste, tendencia, regresión lineal, residuales 
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Introduction 
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Flood risk estimation of a certain area and hydrological dimensioning of 
water works: reservoirs, bridges, walls and protective embankments and 
channelings, is carried out based on so-called Design Floods, which are 
extreme flows of the river associated with low Exceedance Probability. 
The most reliable estimation of design floods is carried out through Flood 
Frequency Analysis (FFA) of the available record or series of maximum 
annual flows. FFA assumes that the random process that generates such 
floods is stationary and therefore ¡ts statistical properties do not change 
over time (Khaliq, Ouarda, Ondo, Gachon, 8: Bobée, 2006; El-Adlouni, 
Ouarda, Zhang, Roy, € Bobée, 2007; El-Adlouni 8 Ouarda, 2008). 


All physical changes that occur in drainage basins, such as the 
construction of  reservoirs, urbanization, channel  straightening, 
deforestation, and land use change alter the hydrological processes, 
giving rise to series or records of non-stationary annual maximum flows, 
because they present trends and changes in variability. In addition, the 
impacts of global climate change can exacerbate hydrological process 
alterations, favoring the formation of non-stationary extreme hydrological 
data sets. Since eany 19905, non-stationarity due to urbanization and/or 
climate change was taken into account, correcting the available records, 
applying incremental factors to predictions, or making predictions in the 
rural basin and then adjusting the results for their level of urbanization 
(Jakob, 2013; Prosdocimi, Kjeldsen, € Miller, 2015). 


Currently, the statistical Extreme Value Theory (EVT) has been 


extended to non-stationary conditions, to fit the classical distribution that 
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follow asymptotically the maximum hydrological data series, ¡.e. the 
stationary probabilistic model Generalized Extreme Value (GEVo) with 
three fitting parameters (u,a,k). This is based on the introduction of the 
so-called covariates, and one of the most used is time (t) in years, 
through which the trend observed in data series can be taken ¡nto 
account, adopting the location parameter variable (us). When the trend is 
linear, U:=po0+p1:t and then the four fitting parameter GVE1 model ¡s fitted 
(Jo, 11,0,K). 1f the trend is curve, ut=po0+pi:t+p2t? and the five fitting 
parameter GVE2 model is applied (uo, u1,142,0,k). In the GVE11 model, both 
the location parameter and the scale parameter (a:+=00+01:t) are linear 
functions of t time and then, in addition to the trend, the change over 
time of series vaniability is considered (Khaliq et a/., 2006; El-Adlouni € 
Ouarda, 2008; Aissaoui-Fgqayeh, El-Adlouni, Ouarda, 8 St. Hilaire, 2009; 
Jakob, 2013; Gado € Nguyen 2016). 


Khaliq et al. (2006), and Meylan, Favre and Musy (2012) have 
presented a review of various methods used in FFA, in the non-stationary 
context. Such techniques include the incorporation of trends in 
distribution parameters or at statistical moments, the quantile regression 
method, local likelihood processes, and the Bayesian approach. In the 
probabilistic models of the first approach, the use of covariates aims to 
integrate changes that have occurred in the past directly into the 
techniques of FFA, with the aim of extrapolating them to the future. In 
this context, models that use time as a covariate, case of the one 
presented here, are simple and useful to reproduce the non-stationary 
behavior of a flood record, but the acceptance of their predictions is not 


entirely correct, given that, as indicated by López and Francés (2013), 
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trends may change in the short or long term, due to climate variability or 
the intensification of human activities. 


Due tothe above, the possibility of incorporating climate indices as 
external forcing into the FFA models has been explored over the past two 
decades, assuming linear or non-linear dependence. This approach has 
shown that the use of such covariates leads to models that better describe 
temporal changes occurred in flood records. In fact, a remote connection 
has been found between the observed changes in flood regimes and the 
anomalies in the indices (ENSO: El Niño-Southern Oscillation; PDO: 
Pacific Decadal Oscillation and NAO: North Atlantic Oscillation) describing 
the temporal evolution of low-frequency atmospheric circulation pattems 
(Khaliq et a/., 2006; López-de-la-Cruz 8 Francés, 2014). 


The objective of this work was to expose the theory of fitting 
Generalized Extreme Values (GEV) distribution in non-stationary 
conditions, by generalizing the L—-moment method, proposed by Gado and 
Nguyen (2016) forthe GVE11 model. Three series of extreme hydrological 
data taken from the specialized literature are exposed and processed, and 
their results discussed, highlighting the simplicity and usefulness of the 
L-moment method to obtain predictions within the historical and future 
record, when they are important, due to the series upward trend and/or 


its increasing variability. 


Operational theory 
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Generalized Extreme Values Distribution 


Annual extreme values are rare events occurring in the right extremity or 
tail of the probability distribution function (PDF), which is the probabilistic 
model that defines the behaviorof the hydrological random variable under 
study, such as: floods and maximum values of rainfalls, temperatures, 
winds and sea levels. The extreme values can be predicted based on the 
PDF, such as the maximum of the random variable that corresponds to a 
certain average interval of recurrence or period of return (Tr), whose 


exceedance probability is p = 1/Tr. 


Extreme value theory justifies and states that extreme data follows 
in an asymptotic manner any of the three types of distributions called: 
Gumbel, Fréchet or Weibull (Clarke, 1973; Stedinger, Vogel, 8 Foufoula- 
Georgiou, 1993; Coles, 2001; Khaliq et al., 2006). These three 
probabilistic models can be represented in one, called Generalized 
Extreme Value distribution (GEV), whose application has been widely 
recommended for modeling annual maximum flows (q) and other extreme 
data (Hosking € Wallis, 1997; Papalexiou €: Koutsoyiannis, 2013; Gado € 
Nguyen, 2016; Stedinger, 2017). The PDF of the GEV with a non- 
exceedance probability [F(q)] is: 
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F(q) = exp(- [1 a whenk + 0 (1) 


In the above expression, u, a, and k are the location, scale, and 
shape parameters of the GEV distribution. When k = O the Gumbel 
distribution is obtained, which is a straight line on the Gumbel-—Powell 
probability paper(Chow, 1964), sothe interval of the variable is: —-oo < q 
< co, When k > 0 the distribution is Weibull, which is a curve with 
downward concavity and upper limit, so: -o <q <u+ a/k. Finally, ifk 
< O the distribution is Fréchet, which is also a curve, but with an upward 
concavity and lower border, so: u + a /k< q <o. The expected 
predictions (Qrr) are obtained with the inverse solution of Equation (1): 


Qu =U + tl — [-In(1 —p)]*) whenk % 0 (2) 


Data sample L-moments 


L-moment method is perhaps the simplest one, and it has become one of 
the reliable procedures for estimating the PDF fitting parameters used in 


hydrology. This is because L-moments, which are linear combinations 
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(Hoking € Wallis, 1997) of weighted probability moments (6), are not 
significantly affected by sample outliers. The first three L-moments of a 
sample (/1, /2, /3) and the asymmetry L-ratio (t3) are estimated by the 
unbiased estimator (b,) of the Br, as follows: 


Eb (3) 
l,=2*b, —b, (4) 
l¿=6-b,-—6-b, +b, (5) 
t=EJE (6) 


The unbiased estimatorof the f,- is (Hoking 8 Wallis, 1997): 


_iym  G-DG-D GN 7 


Tr y Sj=r+1 (n-1)(0-2)--(m-r) + 


where r = 0, 1, 2,.... and gj are the data of the available sample or series 
of size n, ordered from minor to major(q, < q, << q,,). 
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Stationary GEV distribution fitting parameters 


With the L-moment method for the GVEo model, the following equations 
are obtained for its three fitting parameters (Stedinger et al., 1993; 
Hosking 8 Wallis, 1997; Rao 8 Hamed, 2000; Stedinger, 2017): 


k =7.8590 - c+ 2.9554 - c? (8) 
being: 

c= a" 0.63093 (9) 
a = — (10) 


- (1-25.1(1+k) 


u=l, -¿1-Td+1] (11) 


The Stirling formula (Davis, 1972) was used to estimate the value 


of Gamma [(w) function: 
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F(w) = o a (1 4 5 SS 4) (12) 


12w  288:w?  51840:w2  2488320:w* 
Non-stationary analysis approach 


It has two inherent assumptions, the first accepts that non-stationanty 
of the hydrological series of extreme annual values is caused by gradual 
changes in land use or by global or regional climate change; generating a 
slight alteration of its statistical parameters. The second assumption 
accepts that the PDF is time-independent, so the GEV distribution, with 
variable fitting parameters over time, is acceptable for modeling extreme 
non-stationary data. 


Due to non-stationarity, there is great statistical variability 
associated withthe sample or series of available extreme annual data and 
itis therefore appropriate to adopt that mean (1) and standard deviation 
(0) are functions of time. The simplest way, is carried out by means of 
the probabilistic model GEV:1 with linear variations, so: ue=po0+p1:t and 
In(0+)=00+01:t. Inthe above expressions, tis the time-in-years covariate 
that encompasses the n-size data series. The use of natural logarithms of 
Cc ensures positive values of the scale parameter (a) in the GEV 


distribution (Gado 8 Nguyen, 2016). The general expressions ofthe mean 
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and standard deviation of the GEVo model, the second one obtained 
through the second central moment, are (Rao € Hamed, 2000): 


p=u+ 211 T(0+k)] (13) 


0?=SITA+ 2-4) 1214) (14) 


GEV;: distribution fitting 


Since the mean and the natural logarithm of the standard deviation are 
linear functions of time, the expression: ue=po+pi:t is first entered in 
Equation (13) to clear the variable location parameter expression over 
time (us), in which the scale parameter (a+) is also variable, and where 
FK1 isthe first function factor of the form parameter (k), such an equation 
is (Gado 8 Nguyen, 2016): 


U; =p +4 tE [1 TA +] = 9 + 1 0t—FK,- 0; (15) 


where: 
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FK, = 1-T(+k) (16) 


From Equation (14), you get the vanmable scale parameter 


expression over time: 


pe PP (17) 


3 T(1+2:)-12 (1+k) 
So: 
a, = FK, 0, = FK, exp(0, + 0, "t) (18) 


being: 


po razo (14k) (19) 


Substituting Equation (18) in the Equation (15), the variable scale 


parameter expression over time is obtained: 
U,¿ = My +4, t— FK,*FK,"exp(o, +0, *t) (20) 
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The magnitudes Ho and 1 of the straight line representingthe linear 
trend of the sample or data series q:are obtained based on the equations 
of the linear regression straight line, as follows. The dependent variable 
(y) is considered to be the maximum annual hydrological data q: and the 
times or years t are the abscissa (1), in this case equal to the ¡-th value 
of ¡. The slope (11) of the regression straight line adjusted by least squares 
of the residuals and the y-intercept (Jo) are obtained with the following 
equations (Campos-Aranda, 2003): 


q; = Mo + Hot (21) 
1 
_ covíqt) ima qpi-qt 
1 var (t) => 8-2 (22) 
My == Mt (23) 


In the equations above, q and t are the arithmetic means of data 
series and time. The linear correlation coefficient (rxy) measures the 
degree of dependence or association between the variables q and t, varies 
from zero to one, indicating with the one the perfect regression, its 


equationis: 
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_ cov(q,t) 
Txy sm y var(t)"var(q) (24) 


being: 


var(q) ==Ef, q? q? (25) 


To estimate the lineartrend in the standard deviation (0+), first the 
residuals (£e:) are obtained with the following equation (Cundenik 8 Burn, 
2003; Gado 8 Nguyen, 2016): 


El = 9H "t (26) 


Mean value is obtained (s<) and absolute deviations (da) of the 


residuals are calculated, with the following expression: 


ed0 = le, — €l (27) 


Natural logarithms are applied to ef” and represented by a straight 


line whose equation is: 


In(eg%) =0, +0, :t (28) 
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The values of 01 and Oo are calculated based on equations (22) and 
(23). The degree of association between In(<f*) and t is estimated with 
Equation (24) and is designated r,, since the y series is the left side of 
Equation (28). The following equations apply to obtain the stationary 
series without trends in mean and standard deviation (q) (Cundenlik 8 
Bum, 2003; Khaliq et a/., 2006; Gado 8 Nguyen, 2016): 


q” =e8€,—0, fe, >8 when 0t>0 (29) 
ql” =e8e,+0, fe, <ge when 0:>0 (30) 
q” =8,+0, fe, >¿8 when 0;<0 (31) 
q” =8€,—0, fe, <é when 0<0 (32) 


In the above expressions, 0: = exp(00+ 01:t). Equations (9) and (8) 
are applied to the series q5” to estimatethe form parameter (k) of the 
stationary GEVo distribution. Finally, Equation (2) is applied to carry out 
the desired predictions, usingthe expressions of the variable location and 
scale parameters over time, equations (20) and (18), respectively. 
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Standard error of fit 


Since the mid-1970s, the standard error of fit (EEA) was formulated as a 
quantitative measure that estimates the descriptive ability of the fitting 
probabilistic model (Meylan et a/., 2012). The EEA allows objective 
comparison between the various PDFs that are tested or adjusted to a 
series or sample of data, since it has the units of the data (q;). Its 
expression is as follows (Kite, 1977): 


n 6.02 
ERA = [Ei aao* (33) 
n—npa 


where n is the number of data in the available series, npa is the number 
of PDF fitting parameters being tested (five for the GVE::), qí are the data 
sorted from lowest to highest and q, are the values estimated with the 
inverse solution of the PDF (Equation (2)), for the non-exceedance 
probability P(X < x) estimated with the Weibull formula (Benson, 1962): 


P(X <x)= E (34) 


+1 
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in which miis the data order number, with 1 for the lowest and n for the 
highest. 


Approach to probabilistic analysis 


For the GEV:1 distributionto be applicable to non-stationary record, its 
graph of data values (q;) against time (t¡j) must show linear trend and 
variability or standard deviation that increases or decreases over time. 
The calculation of the standard error of fit (EEA) with Equation (33), allows 
the comparison or matching of the various non-stationary models tested 
in the series being processed. When the EEA values are similar, a non— 
stationary model can be adopted in a subjective manner, for example, 


which leads to the most unfavorable predictions. 


For other non-stationary distributions, there are several such as 
Log-Normal (Vogel, Yaindl, 8 Walter, 2011; Aissaoui-Fgayeh eta/., 2009) 
and the Generalized Logistics (Kim, Nam, Ahn, Kim, € Heo, 2015) and 
Generalized Pareto models, which are susceptible to treatment identical 
to that described in this paper for GEV distribution, changing equations 
(13) and (14) to those corresponding to these models (Rao € Hamed, 
2000). Other covariates can also be used instead of time (t), for example 


some global or regional climate indices (López-de-la-Cruz €: Francés, 
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2014; Álvarez-Olguín 8: Escalante-Sandoval, 2016; Campos-Aranda, 
2018a). 


Based on Equation (2), predictions with return periods (Tr) of 2, 10, 
25, 50, and 100 years are estimated through the registration period, 
applying equations (20) and (18). The first prediction corresponds to the 
median, since its non-exceedance probability (1-p) is 50 % and the 
following four are calculated forthe following values: 0.90, 0.96, 0.98 and 
0.99 respectively. In addition, future predictions are made in one of the 
processed series, in years 2025 and 2050. Extrapolating over 30 years of 
the observed behavior in the historical trend and its variability is 
considered to be extremely risky. 


In this regard, mobile L-moment method suggested by Mudersbach 
and Jensen (2010) and applied by Campos-Aranda (2018b), whose main 
limitation is to require fairly extensive records, is considered much more 
reliable. Also shown in data and prediction graphs the estimates of 
extreme Tr (2 and 100 years) with stationary GEVo model, which are 
horizontal dotted straight lines. 


Processed hydrological series 
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Series 1: With a downward trend and a decrease in 
variability 


This record was exposed by Katz (2013) as an impressive example of the 
reduction in the winter (May to October) maximum daily precipitation 
(MDP) measured at the Manjimup station in the extreme southwest of 
Australia. The record covers 75 annual values (1930 to 2004), with a 
declining trend and greater variability atthe beginning of the record. Their 
approximate magnitudes are shown in Table 1, as they were read from a 
bar graph and are shown in Figure 1. 


Table 1. Maximum annual data to be processed for maximum daily 
precipitation (MDP, in millimeters) and flow (q, m3/s) at the indicated 


stations. 
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Figure 1. Time series of data (MDP) and estimated prediction curves 
with GVE:1 distribution at Manjimup pluviometric station, Australia. 


Series 2: With an upward trend and an increase in 
variability 


This annual maximum flows record was processed by Vogel et al. (2011), 
it has 69 values in the period from 1940 to 2008, and comes from a 


gauging station located onthe Aberjona River just afterthe city of Boston, 
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Massachusetts, USA, so ¡its basin has always been under the increasing 
impact of urban development. Their approximate values are shown in 


Table 1, as they were read in a dispersion diagram and are shown in 


Figure 2. 
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Figure 2. Time series of data (q) and estimated prediction curves with 
the GEV:1 distribution at a gauging station on the Aberjona River, USA 
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Series 3: Segment with an upward trend and an 
increase in variability 


This annual maximum flow record was presented by Gilleland and Katz 
(2011), and Katz (2013), from a basin in Washington, USA whose Mercer 
Creek gauging station has 51 values in the period from 1956 to 2006. In 
the 15-year interval between 1971 and 1985, this record was influenced 
by an increasing and intense urbanization. Because of this, only such a 
lapse should be considered non-stationary. The approximate data of that 
series are given in Table 1, as they were read from a bar graph (Katz, 
2013) and are shown in Figure 3. 
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Figure 3. Time senies of data (q) and estimated prediction curves with 
the GEV:1 distribution at Mercer Creek gauging station, Washington, 
USA. 


Description of results 
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Predictions at Manjimup station 


The application of equations (22) to (25) to the 75 MDP data in Table 1 
led to the following results: o = g 52.8382, ui = -0.2811, rxy = -0.4313, 
Co = 2.2767, 01 = -0.0132, and ry = -0.2701. The stationary series 
qs”according to equations (29) to (32) accepts a GEVo model with shape 
parameter: k = -0.0580. 


Equation (2) using the variable location (Equation (20)) and scale 
(Equation 18) parameters leads to a standard error of fit (EEA) of 15.5 
mm. Table 2 shows some of the predictions within the historical record, 
and in the future its predictions are not important due to the downward 
linear trend and the decrease in vanability. The predictions with GVEo 
stationary model of return periods 2 and 100 years are: 38.7 and 96.6 
mm. Figure 1 shows the data and the downward prediction curves. 


Table 2. Predictions (mm) in the historical period at the Manjimup 


station, based on non-stationary GEV:1 distribution. 


Year Data Return periods in years 
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Predictions on the Aberjona River 


The application of equations (22) to (25) to the 69 annual maximum flows 
(q) in Table 1 led to the following results: o = ¿ = 6.1849, 1 = 0.1902, 
xy = 0.4256, 00 = 0.2563, 01 = 0.0304, and r;, = 0.5496. The stationary 
series q” according to equations (29) to (32) accepts a GEVo model with 
a shape parameter: k = -0.1857; instead, Equation (2) using the variable 
location (Equation (20)) and scale (Equation (18)) parameters leads to a 
standard error of fit (EEA) of 2.39 m3/s. Table 3 shows some of the 
predictions within the historical record and in the future their predictions 


are important due to the upward linear trend and the increase in 
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variability. Since this record of 69 data ends in 2008, then the value of 
time t in 2025 is 86 and in 2050 is 111. The predictions with GVEo 
stationary model of return periods 2 and 100 years are: 10.1 and 53.3 
m3/s. Figure 2 shows the data and the upward prediction curves. 


Table 3. Predictions (m3/s) in the historical period and in the future at 


Aberjona River, based on the non-stationary GEV:: distribution. 


No. | Year Data Return periods in years 
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Predictions at Mercer Creek Station 


The application of equations (22) to (25) to the 15 data from 1971 to 
1985, ¡.e., data with numbers 16 to 30 in Table 1, led to the following 
results: Jo = g€ = 6.9048, 1 = 0.5311, hy = 0.6299, 00 = -0.7769, 01 = 
0.1292 y ry = 0.3835. The stationary series q” according to equations 
(29) to (32) accepts a GEVo model with a shape parameter: k = -0.1451; 
instead, Equation (2) usingthe variable location (Equation (20)) and scale 
(Equation (18)) parameters leads to a standard error of fit (EEA) of 0.85 
n?/s. 

Table 4 shows a portion of the predictions within the historical 
record and Figure 3 ¡llustrates the data, the upward curves of predictions, 
and their respective extrapolations. 


Table 4. Predictions (m3/s) in the historical period at the Mercer Creek 
station, based on the non-stationary GEV:1 distribution. 


Year Data Return periods in years 
E mid a 


[2] 


er 
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For this data series, the predictions with the stationary GEVo model 


are of capital importance to approve the extrapolation, toward the 
beginning and end of the record, of the predictions of GEV11 model that 
will be the most critical one (Gilleland 8 Katz, 2011). Therefore, two 
fittings were made: (1) withthe initialstationary record of 15 data (1956- 
1970) and (2) with the lapse of 21 values (1986-2006) occurring after 
the trend and increasing variability segment (1971-1985). The 
predictions of the first fitting were: 5.7, 7.0, 7.6, 7.9, and 8.2 m/s for 
the return periods (Tr) of 2, 10, 25, 50, and 100 years, with an EEA = 
0.193 m/s; the predictions for the second fitting were: 12.2, 19.9, 23.4, 
25.9 and 28.2 m?/s, with an EEA = 0.691 m/s. 


Then, according to Table 4 results, all predictions of GEV:: are 
extrapolable toward the start of the record, being greater than those of 
GEVo. Towards the end of the record, only the median is extrapolated (Tr 
= 2 years), which is 14.2 m3/s, due to the predictions of the other Tr; 
those obtained with GEVo are more extreme than those in the last line of 
Table 4. The above is shown in Figure 3. 
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Discussion of results 


From a practical point of view, the three numerical applications described 
cover the common cases in which is applicable the GVE11 model to non-— 
stationary extreme hydrological data records; these are: 


1. Series with downward trend and variability, present in records of floods 
when several small reservoirs have been built within the drainage basin; 
such behavior may also be associated with global climate change, which 
was the case with the Manjimup pluviometric station. 


2. Records of flows with upward trend, usually associated with basins that 
have intense urban devel-opment or deforestation, which is the case with 
Mercer Creek station. 


If both processes occur inthe basin, or are enhanced by the impacts 
of climate change, variability also increases toward the end of the sample, 


which was the case with the Aberjona River. 


Predictions,beyondthe historical record (t > n), associated with low 
exceedance probabilities (Design Floods) lack importance in records with 
downward trend and variability, as they are smallerin the future. The 
opposite is true for series with upward trend and increasing variability, 
but in these series it is necessary to explain or account for the probable 
physical origin of such behavior in order to accept future predictions and 


try to discern their actual extent. This, because it is extremely risky to 
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accept the behavior of increasing vanability, when, for example, it is not 


known whether urban development will continue in the drainage basin. 


To further explainthe above, itis indicated that Meylan et a/. (2012) 
have noted that the various approaches or new methods of the FFA, which 
accept nonstationarity, have generated techniques that allow predictions 
to be made, associated with low exceedance probability to face the design 
of hydraulic works. However, there is an urgent need to replace the 
concept of retum periods (Tr, in years), to adapt it to the context of 
nonstationarity. This is because the Tr is defined with the average 
recurrence interval, measured over a large number of occurrences. As it 
happens now, in nonstationarity, mean changes over time, the concept of 
Tr is meaningless. Details and solutions have been proposed in the 
references of Sivapalan and Samuel (2009); Salas and Obeysekera 
(2014); Serinaldi (2015), and Salas, Obeysekera and Vogel (2018). 


Conclusions 


Flood Frequency Analysis (FFA) in non-stationary records showing non- 
constant trend and/or variability will be increasingly common in the 
immediate future, due to increased demand for drinking water and food, 


as well as to the impacts of climate change and urban development. A 
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simple and uncomplicated computational approach that allows the FFA to 
be performed on such records is based on the extension of the extreme 
values theory through the fitting with L—-moments, of non-stationary 
distribution GEV11 with location parameter and scale variables lineany 
with the time that is entered as covariate. Therefore, GEV11 model is 
suitable when the series of extreme hydrological data shows trend and 
increase or decrease in variability that can be accepted as linear. 


Based on the results of the three numerical applications in non-— 
stationary records, the simplicity ofthe above method is observed, as well 
as the ease of obtaining the predictions associated with non-exceedance 
probabilities. Graphical contrast is basic for validating the descriptive 
ability of predictions within the historical record and in the near future. 
The numerical results of the standard error of fit will allow the contrast 
and acceptance of other probabilistic non-stationary models. The 
application of the GEV:: probabilistic model to other extreme hydrological 
data, such as winds, temperatures, droughts and sea levels, is fully 
feasible. 
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